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Dissipative dark matter, where dark matter particles interact with a massless (or 
very light) boson, is studied. Such dark matter can arise in simple hidden sector 
gauge models, including those featuring an unbroken U(l)' gauge symmetry, leading 
to a dark photon. Previous work has shown that such models can not only explain 
the LSS and CMB, but potentially also dark matter phenomena on small scales, such 
as the inferred cored structure of dark matter halos. In this picture, dark matter 
halos of disk galaxies not only cool via dissipative interactions but are also heated 
via ordinary supernovae (facilitated by an assumed photon - dark photon kinetic 
mixing interaction). This interaction between the dark matter halo and ordinary 
baryons, a very special feature of these types of models, plays a critical role in 
governing the physical properties of the dark matter halo. Here, we further study 
the implications of this type of dissipative dark matter for disk galaxies. Building 
on earlier work, we develop a simple formalism which aims to describe the effects of 
dissipative dark matter in a fairly model independent way. This formalism is then 
applied to generic disk galaxies. We also consider specific examples, including NGC 
1560 and a sample of dwarf galaxies from the LITTLE THINGS survey. We find 
that dissipative dark matter, as developed here, does a fairly good job accounting for 
the rotation curves of the galaxies considered. Not only does dissipative dark matter 
explain the linear rise of the rotational velocity of dwarf galaxies at small radii, but 
it can also explain the observed wiggles in rotation curves which are known to be 
correlated with corresponding features in the disk gas distribution. 
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1 Introduction 


The evidence for the existence of nonbaryonic dark matter in the Universe has 
grown steadily over the years. In particular, dark matter is needed to explain the 
observed large scale structure and Cosmic Microwave Background radiation (CMB). 
The angular power spectrum of the anisotropy of the CMB has been measured with 
impressive precision by the Planck satellite [1], South Pole telescope [2], WMAP 
[3] and many other instruments. Cold dark matter models, both collisionless [4] 
and even those with nontrivial self-interactions, e.g. [5, 6, 7], provide a simple 
explanation for this data. On much smaller scales, dark matter is also implied by 
measurements of the rotation curves of disk galaxies [8]. Typically such curves are 
asymptotically flat out to the visible edge of galaxies, a behavior in sharp contrast 
to the Keplerian decline anticipated from Newton’s law of gravity. 

Although collisionless cold dark matter does an excellent job in explaining the 
large scale structure and CMB, it is much less successful in explaining the measured 
galactic rotation curves, e.g. [9, 10, 11, 12]. In particular, rotation curve data 
prefer dark matter with a cored distribution rather than the cuspy profile predicted 
from collisionless cold dark matter simulations [13] (for a review and more detailed 
bibliography see also [14]). However, dark matter with nontrivial particle physics 
properties can also explain large scale structure and CMB, and these nontrivial 
properties might play important roles on galactic and subgalactic scales. If this is the 
case, then rotation curves might contain important clues (such as the aforementioned 
preference for cored dark matter profiles) which can help unravel the mystery of dark 
matter. 

It has been argued in previous work [15, 16, 17, 18, 19, 20] that certain dissi¬ 
pative dark matter models can in fact lead to dark matter halos consistent with 
observations. The general idea is that the dark matter halo of a disk galaxy takes 
the form of a dark plasma, with dark matter particles self interacting via a mass¬ 
less dark photon. Such a scenario can naturally arise if dark matter results from a 
hidden sector of particles and forces, where the hidden sector contains an unbroken 
U(l) f (dark electromagnetism) gauge interaction. In such a picture the halo can be 
modeled as a fluid, governed by Euler’s equations. The halo has substantial cooling 
from processes such as dark bremsstrahlung. Heating is also possible. Within this 
framework it was argued that ordinary core collapse supernovae could supply a sub¬ 
stantial halo heat source. This mechanism requires only a tiny interaction coupling 
ordinary particles (electrons and positrons produced in the supernova core) to light 
MeV scale dark matter particles, with kinetic mixing of the photon with the dark 
photon identified as a theoretically favored candidate. This energy is further pro¬ 
cessed into dark photons in the region around the supernova; the end result is that 
a significant fraction of the supernova’s core collapse energy can be transmitted to 
the halo via these dark photons. 

In this scenario the dark matter halo is a complicated dynamical system. Consid¬ 
erable simplification is possible if the system evolves to a steady-state configuration 
where heating and cooling rates approximately balance. It was argued in previ¬ 
ous work that such a picture can explain why dark matter forms a cored profile in 
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disk galaxies. Schematically, the cooling rate from dark bremsstrahlung and other 
dissipative interactions at a particular location, P, is proportional to the square 
of the dark matter number density, n(r), while the heating rate is proportional to 
the product of the number density and dark photon energy flux at that location, 
F 1d (r). The particular process considered being dark photoionization with cross- 
section, < tdp■ That is, 

dT cool = An(r) 2 dV , 

dr heat = n(r)F 7D (r)a DP dV (1) 


where A is a constant which depends on the parameters of the particular dissipative 
model. Given that the dark photon flux produced from ordinary supernovae will 
have some frequency spectrum, the above cross section represents an appropriately 
weighted frequency average. 

Matching the heating and cooling rates of Eq.(l) implies that 


n(r) 


F 1d (y) a dp 
A 


( 2 ) 


At distances r 3> rp (where rp is the exponential disk scale length associated with 
the stars: £* = ^ypre -r ' /rD ) the flux of dark photons falls at the ‘geometric’ rate 

F 1d oc 1/r 2 , while for r ~ rp, F lD falls much more modestly. 2 In fact, simple 
calculations [20] suggest that F lD oc log(r ) for r ~ r D . These features motivate a 
quasi-isothermal dark matter distribution, with dark matter core radius, ro, scaling 
with rp, r 0 ~ rp. 

The dark matter density profile, Eq.(2), scales with the dark photon energy flux. 
This also suggests a correlation with the dark photon source, the type II supernovae 
and by extension, the local star formation rate. This logic suggests that the dark 
matter properties, in particular the density profile, should be closely linked with the 
baryonic properties. The ro ~ rp scaling discussed above is one example of this. 
Another example is a correlation between the supernova rate and asymptotic value of 
the rotation velocity, which has been shown [20] to be equivalent to the Tully-Fisher 
relation [21], Eq.(2) also holds locally: Regions of a galaxy with low gas density can 
have reduced supernova rates leading to a locally depressed dark photon energy flux, 
F 1d . Halo dynamics should also adjust the cooling rate in a corresponding fashion, 
thereby necessitating a lower dark matter density in that region. Thus galaxies with 
wiggles and other features in their baryonic gas distribution can have corresponding 
features in the dark matter distribution, and hence leave an imprint on the rotation 
curve. The purpose of this paper is to explore this possibility. To do so, we utilize 

2 This behavior assumes that the halo is (at least approximately) optically thin. That is, the 
optical depth, r, along the dark photon’s path between its production point and absorption point 
satisfies r ~ 1. In general r depends on the photon’s frequency (energy) so that in practice the 
halo might only be optically thin for a range of dark photon frequencies. Actually this detail may 
not be so important as the energy transported to the halo is dominated by dark photons with 
optical depths r ~ 1. This conclusion is also supported by numerical work [17, 18] which includes 
the effects of finite optical depth and allows for a range of possible shapes for the dark photon 
frequency spectrum generated by supernovae. 
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a very simple formula for the dark matter density distribution, suggested by the 
balancing of the dark plasma cooling and heating processes, described above. We 
then study the implied rotation curves and examine specific examples. 

The outline of this paper, then, is as follows. In section 2 we give some of the 
particle physics details of a candidate dissipative dark matter model. In section 3, we 
develop the necessary tools to describe the dark matter density in terms of baryonic 
properties in this dissipative dark matter approach. In section 4 we examine some 
general features of the predicted dark matter density and rotation curves for generic 
disk galaxies. In section 5 we apply this formalism to some specific examples and 
in section 6 we conclude. 


2 Dissipative dark matter models 


The idea that dark matter originates from a hidden (or dark) sector weakly cou¬ 
pled to the standard model has become quite widely discussed in recent years, e.g. 
[22, 23, 24, 25, 26, 27, 28, 29, 30] (see also [31] for a more detailed bibliography). 
Such models allow for a very rich dark matter phenomenology, especially if the 
hidden sector contains unbroken gauge interactions. If the hidden sector features 
stable matter particles charged under an unbroken U{1)' gauge symmetry, then dark 
matter properties can become extremely interesting. In particular dissipative self 
interactions causing dark matter to cool via the emission of massless dark photons 
become possible. Although the dark sector can be quite strongly self interacting, 
dark sector interactions with electrons, protons, and other ordinary sector particles, 
are relatively weak. In the model considered here, only gravity and a very small 
photon - dark photon kinetic mixing interaction (to be defined in a moment) are 
presumed. 

A very generic dissipative dark matter model was considered in [20] and shown to 
provide a consistent explanation for the CMB anisotropy spectrum and large scale 
structure. In addition it was argued that the same model could also potentially 
explain small scale structure: the cored dark matter distribution in galaxies, Tully- 
Fisher and other scaling relations etc. We refer the reader to that reference for these 
details, here we just review a few essential elements of that model to illustrate the 
basic ideas. 

Consider then a hidden sector featuring two massive stable particles, F l and F 2 . 
These dark matter particles can be either bosonic or fermionic; in the following we 
will assume they are spin 1/2 fermions for definiteness. These particles are presumed 
to be charged under an unbroken U[l)' gauge symmetry, that is they interact with 
a massless 1/(1)' gauge boson - the dark photon. The particle interactions can then 
be described by the Lagrangian: 

^SM T F-d ar k T Fint (3) 

where C$m is the standard model Lagrangian, and Cdark describes the interactions 


3 



( 4 ) 


of the dark sector: 

C dark = —F'^F^ + F^iD,Y - m Fl )F 1 + F 2 (iD- m F2 )F 2 . 

Here, F' = d lt A' u — d u A' is the field-strength tensor associated with the U(l)' 
gauge interaction, A' being the relevant gauge held. Also, the covariant derivative 
is D )j, = + ig Q'A' , where g is the coupling constant relevant to this gauge 

interaction. The dark sector may have much more complexity, additional gauge 
interactions, additional particles etc, but only the low energy effective theory will 
be important for the astrophysical implications discussed in this paper. This means 
that the QED-like Lagrangian, Eq.(4), can actually represent a wide class of hidden 
sector theories. 

In such a theory it is natural to expect dark matter in the Universe to be the 
result of a particle - antiparticle asymmetry, as the symmetric component would 
be expected to efficiently annihilate in the early Universe (at least for much of the 
parameter space of interest). That is, dark matter in the Universe is composed 
of F\ and F 2 particles (with negligible proportion of antiparticles). Dark matter 
asymmetry and local neutrality of the Universe then imply: 

n Fl Q' F i + n F2 Q' F2 = 0 (5) 

where n Fl and np 2 are the number densities of F\ and F 2 respectively. Evidently 
the U(l)' charges of F\ and F 2 , Q F and Q F , must be opposite in sign, but not 
necessarily equal in magnitude. This is, of course, quite analogous to the situation 
with ordinary matter (iq ~ electron, F 2 ~ proton). 

The C, int part in Eq.(3) has yet to be defined. Constraints from gauge invariance 
and renormalizability strongly constrain possible interactions connecting ordinary 
and hidden sector particles. In this particular model, these constraints limit £, nt to 
just the kinetic mixing interaction [32], which leads to photon - dark photon kinetic 
mixing: 

An* = . (6) 

This kinetic mixing interaction effectively provides the dark fermions, iq and F 2 , 
with a tiny ordinary electric charge proportional to the dimensionless parameter e 
[33], 

A tiny kinetic mixing interaction allows ordinary type II supernovae to become 
a huge energy source for the dark sector. Recall that in standard theory nearly all 
of a supernova’s core collapse energy, typically around few x 10 53 ergs, is converted 
into neutrinos in the hot core on a time scale of about 10 seconds. The dozen 
or so neutrino events detected following SN 1987A supports this basic picture [34, 
35]. However there are still substantial uncertainties which means that the energy 
carried away by exotic particles can be up to around 50% of the total core collapse 
energy [36]. It turns out that for kinetic mixing strength e ~ 10~ 9 around half of 
a supernova’s core collapse energy (~ 10 53 ergs) can be converted into light hidden 
sector particles, say iq, if their mass is less than the supernova core temperature, ~ 
30 MeV [36]. 3 [The F 2 dark matter particle can be much heavier, with GeV-TeV 

3 Natural units with fi = c = ks — 1 are assumed. 
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scale preferred [20].] The relevant particle processes include plasmon decay where 
the effectively massive photon propagating in the plasma decays into FiFi pairs. 
Further processes, such as F\F X —> FiFi —> FiFiy D etc, convert this energy 

into dark photons in the region around the supernova, which ultimately escape and 
provide the heat source for the dark halo. This huge heat source, averaging to around 
~ 10 44 ergs per second for the Milky Way (for e ~ 10 -9 ), can potentially compensate 
for the energy lost due to dissipative interactions. In fact, this imbues the dark halo 
with very rich and nontrival dynamics. The halo can contract or expand, and it 
is anticipated that the halo evolves until a steady-state configuration is reached 
where the energy going into the halo sourced from supernovae, as described above, 
approximately balances the energy dissipated. This is the key assumption going into 
the dissipative dark matter scenario advocated here. 

The above particle physics model provides a simple candidate for dissipative 
dark matter. Importantly, the analysis of the present paper is, in many respects, 
model independent - applicable it is hoped to many possible dissipative dark matter 
models. Among these is the mirror dark matter model, where the hidden sector is an 
exact duplicate of the standard model. That allows an exact, unbroken Z -2 discrete 
symmetry to exist which swaps each ordinary matter with a mass degenerate ‘mirror’ 
partner [37]. The matter particles in this case are the mirror electron of mass 0.511 
MeV and mirror nuclei, H', He', O’, Fe ',... and the role of the dark photon is played 
by the mirror photon [31, 23]. Dissipative dark matter models with light rather than 
massless dark photons, or even with light scalar particle(s), might also be possible. 

3 Formalism 


The physical picture for disk galaxies is that they are embedded in a galactic halo 
that takes the form of a plasma 4 composed of dark matter particles [15]. 5 This 
theme has been developed over the last few years in the context of mirror dark 
matter [16, 17, 18, 19] and the simple dissipative model of the previous section [20]. 
Here we build on this recent work, with the aim of describing such dissipative dark 
matter in a fairly model independent manner. 

In the context of the specific hidden sector model described in the previous 
section, the dark matter plasma dissipates energy via various processes which lead 
to dark photon emission, such as dark bremsstrahlung. The halo is also subject to 
heating, sourced by kinetic mixing induced processes in the core of type II supernovae 
in the disk. These processes (7 —> F\ F \, ee —> F\ F \, ee —> FiFi'yD etc) convert up to 
around half of a supernova’s core collapse energy into light hidden sector particles. 
Ultimately this energy is transformed into dark photons in the region around the 

4 There may also be a subcomponent of dark stars. However, the total mass (and distribution) 
of that subcomponent is uncertain and assumed negligible for the purposes of the present analysis. 

5 Elliptical and dwarf spheroidal galaxies form quite a different class, as these galaxy types are 
largely devoid of baryonic gas and show little current star formation activity. The curious reader 
is referred to [20] for a discussion on how these galaxies might fit into this dissipative dark matter 
framework. 


5 



supernova. These dark photons then transport this energy to the halo, where it can 
be absorbed via dark photoionization. 

Using spherical coordinates and setting the presumed thin disk at 9 = tt/2, the 
energy flux of dark photons at a point P = (r, 9, (ft) within an optically thin halo is 
given by: 


F, d (r, 9, (ft) = j dxft 


d r r - 

47r[r 2 + r 2 


S(r,0) 

2 rr sin 9 cos((ft — (ft )] 


( 7 ) 


Here, E(r, (ft) is the luminosity of dark photons per unit area produced on the disk, 
assumed to originate from ordinary type If supernovae. Since supernovae are discrete 
events, this quantity represents an average over a reasonable time period (say last 
10 million years). 

The flux of dark photons will heat the dark halo. The halo also cools via dark 
bremsstrahlung and other processes. The key assumption is that the halo dynami¬ 
cally adjusts (contracts or expands) so that the heating and cooling rates balance at 
any particular point, P. Given that the heating rate is proportional to the product 
of the dark matter density and dark photon energy flux and the cooling rate is pro¬ 
portional to this density squared it follows that p oc F lD [Eq.(l) and Eq.(2)]. That 
is, the dark matter mass density has the form: 


p(r, 6, (ft) = A d<ft dr r 


-‘SN 


>, 4 >) 


47 T[r 2 + r 2 — 2 rr sin 9 cos(<ft — (ft)\ 


( 8 ) 


where A is a constant which depends on the cross section, supernovae dark photon 
energy spectrum etc. In the integrand above, Esn is the type II supernova rate in 
the disk (the proportionality constant in replacing the dark photon luminosity in the 
integrand of Eq.(7) with Esn is also absorbed into A). In principle A also depends 
on the halo temperature at the point P. However, we will assume, for simplicity, an 
isothermal halo so that the constant A is independent of halo location, i.e. r, 9 , (ft. 
[Of course, A may well vary somewhat between different galaxies.] 

The above equation gives a very simple description of the dark matter distri¬ 
bution in galaxies. Furthermore, it is intimately connected with the distribution 
of ordinary baryons through the type II supernova rate per unit area of the disk, 
E,s'Ar(7, (ft). The previous studies [17, 18, 20] assumed that Esjv(r, (ft) could be ap¬ 
proximated by the stellar distribution, taken to be a Freeman disk [38]. This could 
only be a fairly crude approximation as the Freeman disk is a kind of average over 
the past history of the galaxy; the quantity we require is the supernovae rate at the 
present epoch. Since supernovae are the final evolutionary stages of large stars which 
have a short lifespan (typically less than 10 million years) the rate SsAr(r, 0) should 
be correlated with the current star formation rate in the galaxy under consideration. 

Under the simplifying assumption of an azimuthally symmetric disk, the dark 
matter density, p [Eq.(8)], is also azimuthally symmetric. The contribution to the 
gravitational acceleration at a point in the plane of the disk due to the dark matter 
distribution, Eq.(8), therefore depends only on the distance r and is given by 


g(r) = Gn d(ft dcos9 / dr r 2 


p(r, 9) cosca 


r 2 + r 2 — 2 rr sin 9 cos (ft 


( 9 ) 
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Here Gn is Newton’s constant and 

r — r sin 6 cos 0 

costa = _ . (10) 

y r 2 + r 2 — 2rr sin 6 cos 0 

In this azimuthally symmetric case, the motion of the baryons in the disk can be 
assumed to be circular and the dark matter halo contribution to the rotational 
velocity is then given by v\ aio /r = g(r). 

If the halo does indeed evolve to a steady-state configuration where heating and 
cooling rates locally balance then the distribution of the dark matter out to some 
distance, Thaio, should be fixed by Eq.(8). This means that for 0 < r ~ Thaio a 
kind of equilibrium configuration ensues so that the dark matter distribution is 
approximately independent of the initial conditions and early history of the galaxy. 
6 In particular, the distribution is also independent of the total dark matter mass 
of the halo, m c i a rk■ However, the size of fhaio is expected to depend on this quantity: 

m dark = J pdV ~ ^p^r ha i 0 (11) 

where we have replaced p with the quasi-isothermal distribution, an approximation 
roughly valid, especially at large distances, r > ro (as will be seen in section 4). 
Observations reveal that rotation curves are generally flat at the largest observed 
distances in a given galaxy, which indicates that r ka io is typically greater than this 
observed edge of the galaxy. 


4 Generic Galaxies 


We now examine the case of a ‘generic’ disk galaxy. Specifically we consider a thin 
disk galaxy (with no bulge) with stellar mass distribution described by a Freeman 
disk [38]: 


S*(r) = 


rrio _ t— 

- - Tre r D 


where uid is the stellar mass of the disk and To is the disk scale length. To solve 
Eq.(8) we require the type II supernova rate in the disk, 0). Here we shall 

assume azimuthal symmetry and take E 5 JV (t) oc E*(r). Of course this can only be 
a very crude approximation, but can nevertheless provide useful insight for generic 
disk galaxies. With this assumption, the type II supernova rate per unit area is 
then: 


yW-rl 

£ Sn(t ) = — Rsn 

m D 


6 Regarding the early history of the halo evolution, we expect that at very early times (t ~ Gyr) 
the dark matter that seeded the galactic halo would have collapsed, possibly into a disk. Shortly 
thereafter, following the collapse of ordinary matter and consequent star formation activity, the 
dark matter was reheated. This reheating expands the dark matter gas component into the plasma 
halo and eventually the system should relax into a steady-state configuration where the heating 
and cooling rates locally balance. This is the presumed state of disk galaxies at the current epoch. 
See also [31] for some more detailed discussion. 
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where Rsn is the total supernova rate of the generic galaxy under consideration. We 
further assume that this total supernova rate scales with stellar mass via: Rsn oc 
y/mD, which is consistent with the supernova study [39]. 

With these baryonic relations, it is straightforward to numerically solve Eqs.(8,9) 
to obtain the dark halo density distribution and rotation curve. In doing so, we 
assume that the coefficient A is a constant, which we set by fixing the maximum 
rotational velocity of a = 10 n m© generic galaxy to be 250 krn/s (a typical value 
for such a galaxy). 7 

Consider first the resulting dark matter distribution, Eq.(8). This distribution 
is not spherically symmetric, due to the spherically asymmetric heating from the 
disk. In figure 1 we plot p(r,9) for 6 = 0 (direction normal to the disk), 6 = 7t/6, 
6 = 7t/4, 6 = 7t/3 and 6 = 27t/ 5. The density has a log divergence as 6 —» 7t/2 
(plane of the disk). As the figure shows, the dark matter halo approaches spherical 
symmetry for r 3 > Td, with significant departures from spherical symmetry in the 
inner region (r ~ ryj). 



r(kpc) 

Fig. I: Dark matter density, p(r,9 ) derived from dissipative dynamics [Eqs.(8,13)] for 
various 9 values. Curves from bottom to top are for: 8 = 0 (direction normal to the disk), 
6 = 7t/6, 6 = 7t/4, 9 = 7t/ 3 and 9 = 27t/ 5. The reference values rno = lO 11 ???.©, ryj = 4.3 
kpc are assumed. 

7 In realistic models we expect some variation of A between galaxies. In fact, A oc 1 /y/Thaio 
follows in models where dark bremsstrahlung is the dominant cooling process. Simple analytic 
calculations suggest that the halo temperature is related to the asymptotic rotational velocity via 
Thaio ~ fhitf , ot /2 where m is the mean mass of the dark matter particles forming the plasma [16]. 
Utilizing the Tully-Fisher relation we then have A oc (L^ffi 1 / 4 , where Lb is the B-band luminosity 
of the galaxy in question. 





Fig. 2: Generic rotation curves predicted from this dynamics, Eqs.(8,9,13). The four 
curves are for parameters: (a) rno = lO 8 m 0 , ro = 1.0 kpc, (b) mo = lO 9 m 0 , ro = 1.4 
kpc, (c) rrio = 10 lo m 0 ,rc> = 2.1 kpc, and (d) mo = 10 11 m 0 , ro = 4.3 kpc. The curves 
represent only the dark matter contribution to the rotational velocity. 


We now turn to the rotation curves, Eq.(9). Figure 2 shows the result of nu¬ 
merically solving this equation for four illustrative examples with stellar masses of 
mo = lO 8 m 0 , mo = 10 9 m & , mo = lO lo m 0 and rrio = 10 n m Q . The correspond¬ 
ing values for the disk scale length (given in the figure caption) are typical of high 
surface brightness disk galaxies. 

In figure 3 we compare the predicted rotation curves (dark matter contribution 
only) for two galaxies with the same baryonic mass mo = lO lo m 0 but with two 
different disk scale lengths: (a) ro = 2 kpc and (b) ro = 7 kpc. These parameters 
correspond roughly to the high surface brightness galaxy NGC 2403 (a) and low 
surface brightness galaxy UGC 128 (b). [To do an actual fit to the rotation curve 
data for these two galaxies, a more sophisticated treatment of T,sn is recommended, 
such as the approach in section 5 which makes use of a Kennicutt-Schmidt-type 
relation.] Although both galaxies have the same asymptotic value for the rotational 
velocity, the low surface brightness galaxy has a much shallower rotation curve. This 
behavior has been known for many years, see in particular the discussion in [40, 41], 
and is very simply explained with dissipative dark matter dynamics. 
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0 10 20 30 40 50 

r (kpc) 

Fig. 3a: Predicted rotation curve (dark matter only contribution) for two galaxies with 
the same baryonic mass mp = 10 10 mQ but with two different disk scale lengths: (a) 
ro = 2 kpc and (b) ro = 7 kpc. 



r/r D 


Fig. 3b: The same two curves of figure 3a plotted in terms of the dimensionless variable 
r/rr>- The curves coincide as rr> is the only length scale in the problem. 
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Phenomenologically, it has long been observed that spherically symmetric cored 
dark matter distributions provide a reasonable fit to the rotation curves of typical 
disk galaxies, e.g. [42, 43, 44, 45]. In this context, the spherically symmetric quasi- 
isothermal distribution: 


Piso( r ) = ^ (14) 

r r o 

with rotation velocity 

— = I# f piso (f) 4 nf 2 df (15) 

r r Jo 

is often considered. Other closely related cored profiles have also been discussed in 
the literature, such as the Burkert profile [46, 47]. 

Although the dark matter density profile resulting from the assumed dissipative 
dynamics is not spherically symmetric, it turns out that the dark halo contribution 
to the rotational velocity can be approximated by the above spherically symmet¬ 
ric quasi-isothermal distribution. We illustrate this in figure 4, which compares 
the generic rotation curves of figure 2 [derived from Eq.(8) and Eq.(13)] with the 
rotation curves resulting from the spherically symmetric quasi-isothermal distribu¬ 
tion, Eq.(15). Furthermore, the approximate agreement of the two curves requires 
r 0 ~ r D and p 0 ~ 7 x 10~ 2 m Q /pc 3 , independently of mr> over the considered range: 
10 8 < rriD/iTiQ < 10 11 . This value for ro is somewhat below best fit values, which 
prefer ro ~ 2 m [48]. We suspect that this is due, at least in part, to our use of E* 
to model E sn- In the next section we will use the Kennicutt-Schmidt-type relation 
which we find does a reasonable job in explaining rotation curves. 

Also, the result that po was found to be roughly constant relied on A being a 
galaxy independent quantity and also the assumed Rsn oc yJmE scaling. In models 
where dark bremsstrahlung dominates the halo cooling, A oc 1 /\f Thaio, where Thaio 
is the temperature of the dark matter halo. Simple analytic calculations as in [20] 
suggest Thaio oc Po' r o- ^ I s taken into account, then we would have po oc 1 /r \/ :i , 
a galactic scaling relation that is roughly consistent with observations [49, 50] (the 
Tully-Fisher relation [21] also follows from similar arguments [20]). These scaling 
relations, though, are limited by uncertainties such as that of the Rsn oc y/rrTR 
scaling. 
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r (kpc) 


Fig. 4: Comparison between the rotation curves from dissipative dynamics for the exam¬ 
ples in figure 2 (sold lines) with that using the quasi-isothermal profile [Eq.(14)] (dashed 
lines). 


5 Specific examples 


Over the last few years, high quality rotation curve data has become available from 
surveys of nearby galaxies, e.g. [51, 52], ffere we shall examine examples contained 
in the LITTLE THINGS survey of 26 nearby dwarf galaxies [52], Additionally, we 
consider also the small spiral galaxy, NGC 1560 [53, 54], The selected sample mostly 
feature flat or rising rotation curves out to the largest measured distance. 

The study [52] compared the rotation curves of the LITTLE THINGS sample 
with the quasi-isothermal profile, Eq.(14) and also the NFW profile [13] 

PNFw(r) = ———2 

(r/r s )(l +r/r s ) 2 

It was found that the NFW profile generally gave very poor fits to the observed 
rotation curves, while the quasi-isothermal profile generally fared much better. As 
discussed in the previous section, the quasi-isothermal profile was historically in¬ 
troduced purely phenomenologically, without reference to any dark matter theory, 
but nevertheless can be viewed as a zeroth order approximation to dissipative dark 
matter (figure 4). Note though, that dissipative dark matter provides an even more 
constrained description since the core scale length (ro) correlates with the disk scale 
length ro, via the association of the heat sources with core collapse supernovae in 
the disk. The purpose of the current paper is to illustrate these ideas with specific 
examples, and postpone a more detailed study to future work [55]. 


(16) 


12 










As in section 4 we shall derive the rotation curves using Eqs.(8,9) for the dark 
matter contribution. However, instead of using £* for the supernova rate in the 
disk, we shall use a Kennicutt-Schmidt-type relation: £sw(r) oc [E 9as (r)] w . Such a 
relation allows E SN to be related to local properties of each galaxy. Explicitly, we 
use 

P( r > 0) = A' J d4> J df f j . -~r (17) 

j j 47 t \r z + r z — Zrr sin t) cos <p\ 

where the coefficient is now A', having absorbed the proportionality constant of the 
Kennicutt-Schmidt relation into its definition. For N = 2, A' has units pc/m & . 

For the exponent, N, in the Kennicutt-Schmidt type relation, we initially consid¬ 
ered the range: 1.0 < N < 2.0, but for the sample of galaxies considered here it was 
found that the rotation curve data could be reasonably well explained if N — 2.0 
(although lower values, e.g. N = 1.5 also gave reasonable results). The value of 
N = 2.0 is consistent with the star formation rate found by Schmidt [56], who con¬ 
sidered star formation rates in the Milky Way, but somewhat higher than the value 
obtained by the more recent study of Kennicutt, N ^ 1.4 [57] (see also the recent 
review [58]). However the Kennicutt study considered the total star formation rate 
in units m & yr -1 averaged over the disk for a given galaxy, while the quantity of 
interest in the present context is the local formation rate of large stars, m star ~ 8m Q , 
the progenitors of type If supernovae. Anyway, with the exponent fixed at N = 2.0, 
we could attempt to fit each galaxy with just one parameter, A'. [Although we allow 
A' to vary between different galaxies, it is found that there is only a modest varia¬ 
tion, within a factor of two or so between the 10 galaxies considered, with central 
value around A' ~ 0.001 pc/m & .] 

We first consider the small spiral galaxy, NGC 1560 [53, 54], This example is 
particularly interesting because it shows a small ‘wiggle’ in the rotation curve, at 
around r ~ 5 kpc. The gas density for this galaxy is given in [54], reproduced in 
figure 5a, also has this feature at r « 5 kpc. Although the wiggle is quite small 
in the azimuthally averaged rotation curve data, it can be seen more easily if one 
separately examines the northern (receding) and southern (approaching) sides of the 
galaxy [53, 54], For this reason, these and other authors, e.g. [41] have argued that 
the wiggle is significant and demands an explanation. 

For the purposes of the present analysis, though, we shall for simplicity assume 
azimuthal symmetry. In solving Eq.(17), we therefore set T, gas to be the average 
of the northern and southern gas density measurements, T, gas (r) = [E^,° rt/l (r) + 
E g °s th ( r )]/2. In figure 5b we compare the predicted rotation curve with the data. 
As the figure shows, dissipative dark matter does a good job in reproducing both 
the rotation curve (Xred = 0.45) as well as the small wiggle at r ~ 5 kpc. The 
phenomenological cored halo models also do a good job in fitting the rotation curve, 
with Xred = 0-30 for quasi-isothermal and Xred = 0-33 for Burkert profile. [The NFW 
profile, on the other hand, gives a much poorer fit to the rotation curve data [54].] 
However such smooth profiles cannot explain the apparent wiggle at around r « 5 
kpc. 

We next consider a selection of dwarf galaxies from the LITTLE THINGS survey 
[52], That reference contains both the gas density and rotation curve data. The 
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results of the analysis are shown in figures 6-14 for the galaxies: DD087, DD0126, 
DD0133, DD052, DD047, DDO70, NGC 3738, NGC 2366 and DD0154. The 
figures show reasonable agreement with the data, especially as only one parameter, 
A', was adjusted. The approach here also seems capable of explaining the subtle 
wiggles in the rotation curves and could easily be improved by taking into account 
the azimuthal dependence of the gas density. The galaxy NGC 1560, where the 
wiggle in both the gas density and rotation curve at r ~ 5 kpc is much more 
pronounced in the northern half of the galaxy [53, 54], might be a good example to 
use if one were to take into account the azimuthal dependence. 

The results in figures 5-14 assumed that the Kennicutt-Schmidt type relation 
had the exponent fixed at IV = 2.0. Allowing N to vary with values greater than 
2.0 generally gave better agreement between the predicted rotation curves and the 
data. We illustrate this in figure 15, where we consider the predicted rotation curve 
for DD0154 with varying N values (where in each case we adjusted X' to give the 
correct normalization). Given the uncertainties in the (large) star formation rates, 
large N values such as IV = 3.0 might be possible. However, we do not expect 
Eq.(17) to be an exact description of the dark matter density profile. In particular, 
the assumptions of a) azimuthal symmetry b) isothermal halo and c) negligible dark 
stellar contribution to the rotation curve, are all expected to be violated, at least 
to some extent, in specific cases. Given these potential uncertainties the level of 
agreement of the predicted rotation curves with the data is reasonable. 

In addition to the galaxies shown, we also considered DDO101 and DD0168, 
which were found to give much poorer fits and are not shown. DDO101 has a rela¬ 
tively low baryonic gas density and total baryonic mass. The resulting suppressed 
heating of DDOlOl’s halo may have been insufficient to support the halo from col¬ 
lapse. That is, the dark halo of DDO101 is no longer a diffuse plasma, but may have 
collapsed into dark stars. This is also supported by the relatively high measured 
value for the central density of its halo. By contrast, DD0168 has a relatively high 
baryonic gas density in the central region. It is unclear why the formalism might 
fail for DD0168 but it might have something to do with this detail. Also, the rota¬ 
tion curve data for both DD0168 (as with DDO70) show a turn over at the largest 
measured distances, which would require an understanding of the boundary region, 
beyond the scope of this work. 
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Fig. 5a: Baryonic gas density as found in [54]. The blue circles (black squares) represent 
an average over the northern (southern) half of the galaxy. 
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Fig. 5b: Predicted rotation curve for NGC 1560 [solid line]. Also shown are the separate 
contributions from the assumed dissipative dynamics, Eq.(17) [dashed-dotted line], the 
baryonic stellar [dashed line] and baryonic gas contributions [dotted line]. The data is 
from [54], 
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Fig. 6: Predicted rotation curve for DD087 [solid line]. Also shown are the separate 
contributions from the assumed dissipative dynamics, Eq.(17) [dashed-dotted line], the 
baryonic stellar [dashed line] and baryonic gas contributions [dotted line]. The data is 
from [52], 
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Fig. 7: Predicted rotation curve for DD0126 [solid line]. Notation as in Fig. 6, data from 
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Fig. 14: Predicted rotation curve for DD0154 [solid line]. Notation as 
from [52], 
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Fig. 15: Predicted rotation curve (dark matter + baryons) for DD0154. 
Schmidt exponent, N, varying. Curves from (left) bottom to top are ft 
1.5. N = 2.0. N = 2.5. N = 3.0. 








6 Conclusions 


We have considered dissipative dark matter, in which dark matter particles form¬ 
ing galactic halos are presumed to have dissipative self-interactions. Hidden sector 
gauge models featuring dissipative dark matter are straightforward to construct. An 
example being hidden sector dark matter charged under an unbroken 1/(1)' gauge 
symmetry, leading to dark matter interacting with a massless gauge boson - the 
dark photon. The dark photon not only mediates self interactions among the dark 
sector particles, but allows for halo cooling via dissipative particle processes such 
as dark bremsstrahlung. Mirror dark matter, where the hidden sector is an exact 
duplicate of the standard model sector, is a constrained example of such a theory. 
Many other dissipative dark matter models may be possible, e.g. where the dark 
photons are light rather than massless, or even where they are replaced by light 
scalar particle(s). 

If dark matter is dissipative then the dark matter halo around disk galaxies can 
have quite nontrivial behavior. In addition to cooling, the halo can have heating 
powered by ordinary type II supernovae in the disk. Naturally, this requires a 
small coupling between the dark matter particles and ordinary matter ones, with 
the photon - dark photon kinetic mixing interaction a likely suspect. With such an 
interaction ordinary supernovae become a major source of dark photons which can 
supply a substantial amount of energy to the halo. 

This type of dark matter, which can be modeled as a fluid satisfying Euler’s 
equations of fluid dynamics, will expand and contract in response to these heating 
and cooling influences. In the current epoch, the halo is anticipated to have evolved 
into a steady-state configuration where the energy radiated at each point approxi¬ 
mately balances the heating. If this indeed occurs, then the physical properties of 
the halo will be strongly constrained and a very simple formula, Eq. (8), relating the 
dark matter density with the supernova energy sources results. 

For a given supernova distribution, a prediction for the rotation curves of individ¬ 
ual galaxies follows - all of the model dependence, kinetic mixing parameter, relevant 
cross-sections etc can be condensed into a single parameter. The spatial distribution 
of supernovae was modelled with a Freeman disk profile and also via a Kennicutt- 
Schmidt-type power law. The resulting rotation curves were then compared with 
observations and found to be in reasonably good agreement. In particular, subtle 
‘wiggles’ in the rotation curve which appear to be correlated with corresponding 
features in the disk gas distribution could plausibly be explained in this approach. 
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